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ABSTRACT 

o 

Context. Mergers of neutron stars (NS) and black holes (BH) are among the strongest sources of gravitational waves and are potential 
' central engines for short gamma-ray bursts. 

' Aims. We aim to compare the general relativistic (GR) results of other groups with Newtonian calculations of models with equivalent 

parameters. We vary the mass ratio of the NS to BH and the compactness of the NS. The mass of the NS is 1.4 Mq. We compare the 
dynamics in the parameter-space regions where the NS is expected to reach the innermost stable circular orbit (ISCO) before being 
tidally disrupted (mass shedding, MS), and vice versa. 
fT^ ' Methods. The hydrodynamics is evolved by a Newtonian PPM scheme with four levels of nested grids. We use a polytropic EoS 

(r = 2), as adopted in the GR simulations. However, instead of full GR we use a Newtonian potential supplemented by a Paczynski- 
Wiita-Artemova potential for the BH, both disregarding and including rotation of the BH. 
' Results. If the NS is compact (C = 0.18), it is accreted by the BH more quickly, and only a small amount of mass remains outside 

^J^ ' the BH. If the mass ratio is small (2 = 2 or 3) or the NS is less compact (C = 0.16 or less), the NS is tidally torn apart before being 

; I accreted. Although most of the mass is absorbed by the BH, some 0.1 Mq remain in a tidal arm. For small mass ratios (2 = 2 and 

. 3), the tidal arm can wrap around the BH to form a thick disk. When including the effects of either BH spin-up or spin-down by the 

CIh' accreted matter, more mass remains in the surroundings (0.2-0.3 Mq). 

Conclusions. Although details and quantitative results differ, the general trends of our Newtonian calculations are similar to the GR 
calculations. A clear delimiting line separating the ISCO from the MS cases is not found. Inclusion of BH rotation as well as sufficient 
numerical resolution are extremely important. 
C/2 . 
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^ ' 1. Introduction We compare the results of a partially post-Newtonian calcu- 

OO , lation of NS-BH mergers with published GR simulations, to ad- 

■ dress in particular the following aspect: Miller (2005) pointed 

■ Merging binaries of neutron stars (NS) and black holes (BH) are "^^^ that during the NS-BH binary orbital inspiral, the inner- 
^ • of interest both as possible central engines for short gamma-ray "^o^t stable circular orbit (ISCO) can be reached before the 
^ ; bursts and as promising gravitational wave sources detectable by ^S becomes tidally disrupted; this was confirmed by Taniguchi 
S ' the interferometer observatories. As computational power and f (^008) and Ferrari et al (2009). In this case, it is to 
2; numerical schemes have progressed so has the detail of Simula- expected that most of the NS matter is swallowed by the 

tions of these mergers. Eai-ly simulations by our group (Jankaet ^^foi"^ *e material can spread out and form a massive 

> ■ al. 1999) were developed to investigate gamma-ray burst see- t"™^- ^mce the ISCO is a purely genera relativistic eff^ect, its 

■ narios. We included detailed microphysics (equation of state, consequences will not be seen in a purely Newtonian simula- 
<S ■ neutrino source terms) but accounted for general relativistic ef- ^i""' However, the Paczynski- Wiita potential (Paczynski & Wiita 

c3 fects either only partially (gravitational wave emission and back- 1980 Abramowicz 2009) phenomenologically mimics the ISCO 

■ ■ reaction) or phenomenologically (Paczynski- Wiita potential for ^ Newtonian setting. So the question arises as to whether and 

the BH). Other groups interested in gravitational wave aspects ^^at extent one can reproduce the results of the general rela- 

have published results of NS-BH mergers with general relativis- ^i^i^^ic simulations with crude Newtonian approximations, 

tic (GR) physics included to different levels of approximation Which of the two cases of NS reaching ISCO first or NS 

(Etienne et al. 2009; Shibata et al. 2009; Rantsiou et al. 2008; becoming tidally disrupted first happens during the evolution 

Taniguchi et al. 2008; Duez et al. 2008), but with some vari- will have a direct impact upon the amount of mass remaining 

ant of the polytropic equation-of-state (EoS) for the neutron star in the torus around the BH. This mass is supposed to provide 

matter Rantsiou et al. (2008) did not include dynamic GR in that the energy for the jet of the short-duration gamma-ray bursts, so 

they used a fixed background. Ferrari et al. (2009) investigated a lower mass should produce a dimmer burst. Intrinsically, the 

quasi-equilibrium sequences of ellipsoidal NSs also including tidal disruption — called 'mass shedding' (MS) in previous pa- 

microphysical EoSs in the tidal approximation. A good source pers — should yield a more massive torus than the ISCO case, 

of references are the reviews of Faber (2009) and Duez (2009). As shown in our previous investigations (Ruffert et al. 2001 and 
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Table 1. Key initial model parameters and some results. 
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Initial parameters: mass ratio Q, compactness C, black hole mass Mbh, neutron star radius i?Ns, initial orbital distance rf;, size of largest grid 
L, size of finest zone A.t = L/2048. The mass of the NS is 1.4 Mq in all cases. Values at the end of the simulation, only for high-resolution 
models: unbound+ejected gas mass M^, bound neutron star mass around BH Mb, neutron star mass not instantly accreted by BH Mg, |h/l| 
states whether the difference between values for M^ in low-resolution and high-resolution simulations is less than 20% (Yes/No; '-' no low-res 
model computed). The spin-up of the BH by accretion of matter is not taken into account. 



references therein, Setiawan et al. 2006), fairly massive tori of 
at least several hundredths of a solar mass seem to be needed for 
the scenario of a hot neutrino-emitting disk to function as a cen- 
tral engine for gamma-ray bursts (Lee et al. 2005, Oechslin & 
Janka 2006). Thus the question of how the NS becomes tidally 
disrupted or swallowed by the BH, and how much mass remains 
inside the torus is an important one to explain gamma-ray bursts 
that needs to be investigated in detail. The results of previously 
published GR simulations still disagree with each other in this 
point. 

Following the introduction above, this paper will describe 
the methods in the next section and the results of the simula- 
tions in Sect. |3] These results are then discussed (Sect. HJi and 
the paper concludes with a summary (Sect.|5]l. 

2. Methods 

2. 1 . Theory and numerical formulation 

The three-dimensional hydrodynamic simulations were per- 
formed with a basically Newtonian code based on the piece- 
wise paraboUc method (PPM) of Colella & Woodward (1984) 
with four levels of nested grids (Ruffert 1992). The equidistant 
Cartesian grids all have the same number of zones, depending 
on the model, either 128 or 256 zones per dimension. Initially 
the finest grid covers both components of the system (neutron 
star and black hole). We use an ideal gas equation-of-state (EoS, 
P-{r- l)e with adiabatic index F = 2, pressure P, and internal 
energy density e) to close the energy equation evolved in the hy- 
drodynamics code. The self-gravity of the gaseous mass distri- 
bution on the grids are calculated using fast-Fourier transforms 
(FFT). Details of the implementation of the nested grids can be 
found in Ruffert (1992). A low density medium, of roughly 10-^ 
the central density of the NS, fills the volume outside the NS and 
the BH. This "atmosphere" is necessary for numerical reasons 
(e.g., use of the reciprocal density) and is common to grid-based 
codes, both for Newtonian and GR simulations. The density of 
all zones is kept above this minimal value, which is referred to 
as "vacuum". 

The black hole (BH) of mass Mbh is treated as a gravita- 
tional point mass surrounded by a vacuum sphere discretised on 




12 3 4 

r/Rs 

Fig. 1. Gravitational acceleration as a function of distance r, as 
given by the prescription of Eq. ([T]l, for various spin parame- 
ters a = 0.999, 0.99, 0.9, 0.5, 0; soHd Hnes from left to right. 
The dotted line shows a Newtonian potential acceleration ccr^^. 
The bold solid line shows the Paczynski-Wiita case oc (r - ^s)^, 
i.e., a-Q. is the Schwarzschild radius. 

the grids. The existence of the innermost stable circular orbit 
(ISCO) is mimicked by including a Paczynski-Wiita-type po- 
tential (Paczyriski & Wiita 1980, Artemova et al. 1996). The 
Artemova et al. prescription for the acceleration is 

dpBH ^ GMbh 
dr r^-P{r - ruf ' ^ ^ 

where th is the black hole event horizon, and /? is a constant 
for a given value of the BH spin parameter a. It is defined by 
/5(a) = (nn(a) - ruia))/ ruia), where rin(a) is the radius of the 
ISCO. This prescription reduces to the usual Paczynski-Wiita 
potential when the BH spin parameter a is zero. We then have 
that 13-2, and rii(a = 0) = 2GMbh/c^. We note a mathematical 
discontinuity: for y6 > 0, re remains a lower Umit for r > th, 
whereas for y6 = one recovers the Newtonian case for all r > 
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(see Fig. [T]i. The graphs of ru{a), /3(a), and rin(fl) can be found 
in Setiawan et al. (2006), their Fig. 1. 

This potential is implemented numerically as follows. First 
the total Newtonian potential is calculated for all gas on the 
grids, including the BH mass. From this potential, the acceler- 
ations are derived by discrete differencing, as is common proce- 
dure. We then add the acceleration given by the difference be- 
tween the formalism of Eq. ([T]i and a purely Newtonian poten- 
tial. By construction, this difference only acts on the gaseous 
NS matter on the grid. To respect Newton's Third Law ("ac- 
tion=reaction"), a tab of the sum of these accelerations has to 
be kept and its inverse enforced on the BH. 

The radius of the vacuum sphere mimicking the BH is cho- 
sen to be the arithmetic mean of the event horizon and the ISCO. 
The mass, momentum, and angular momentum of material flow- 
ing across this inner boundary are added to the BH values and the 
zones affected are reset to near-vacuum values. The position of 
the BH is adjusted accordingly — the centre-of-mass of the gas 
leaving the grid and of the BH do not necessarily coincide — and 
moved by a leapfrog procedure (alternate position and velocity 
updates). In all models listed in Table[T] the angular momentum 
of the accreted material does not change the spin parameter a in 
the equations that determine the BH potential, Eq. ([TJ; i.e., the 
BH potential only deepens because of the increase in mass of the 
BH. We evolved an additional set of four models to investigate 
the effects of BH spin. These models are listed in Table |2] and 
discussed in separate paragraphs. 

The general relativistic effect that is modelled phenomeno- 
logically by the procedures described above is the presence of 
the ISCO and its change in radius. Frame-dragging and many 
other effects, e.g., relativistic kinematics, redshift, time dila- 
tion, are omitted. However, we do include the local effects of 
gravitational-wave emission — the volume integral of which re- 
produces the quadrupole approximation — and the correspond- 
ing back-reaction on the hydrodynamic flow (for details see 
Ruffert et al. 1996). 

2.2. Initial conditions 

The initial distance of NS to BH varies from model to model, to 
ensure roughly 2-3 orbits before coalescence (see Table [TJ. For 
the finest grid to cover both components of the binary, the geo- 
metric extent of the grids has to be adapted to ensure the highest 
possible numerical resolution. Both the extent of the largest grid 
as well as the size of the finest zone are listed in Table [T] 

Parameters for the NS model are given as follows. The 
mass is kept constant for all polytropic models, Mns - 1.4-Mq. 
The radius, /?ns, is then fixed by the chosen compactness, C - 
(GMns)/(c^^ns)- For a F = 2 polytrope, the relation between 
pressure P and density p is P = Kp^, where k = f G/?^^ is fixed 

by the NS radius. The density distribution is p(^) - pc^^^, with 

Mm = t^^lsPc, and ^ = nr/R^. 

The main variable of the BH, its mass, is defined in terms of 
its mass ratio with respect to the NS, Q = Mbh/A^ns- Both the 
NS and BH are initially given Keplerian velocities plus a small 
radial velocity component from the quadrupole approximation 
of gravitational wave emission. 

In our Newtonian framework, there is no ambiguity about 
which masses and radii to use in the equations above. This dif- 
fers from relativistic calculations, where a choice of masses (e.g., 
ADM, rest), and distances (e.g., isotropic, circumferential) ex- 
ists. We return to this point when trying to compare our results 
with other people's simulations, in Sect.|4] 
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Fig. 2. Parameter space of models presented in this paper (only 
non-spinning BH), with compactness C - (GMns)/(c^^ns) and 
mass ratio Q = Mbh/A^ns- Each symbol represents one model: 
'diamond' Mg < IO^^Mq, 'triangle' IQ-^Mq < Mg < O.O2M0, 
'square' O.O2M0<Mg<O.lM0, 'star' O.IMq < Mg, for the mass 
Mg not instantly accreted by the BH. Note that the spin-up of 
the BH by the accretion of matter is not taken into account. 
The line separating the innermost stable circular orbit (ISCO) 
and the mass shedding (MS) regions is taken from Taniguchi et 
al. (2008). 
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Fig. 3. Results of Shibata et al. (2009) and Etienne et al. (2009) 
(only non-spinning BH) models, with compactness C and mass 
ratio Q. Each symbol represents one model; we use the same 
symbols as used in Fig. |2] The symbols for the Etienne 
et al. (2009) models are sHghtly shifted to C = 0.147 (and 
coloured), for clarity. 



When the NS is evolved in the simulation, the initially spher- 
ical NS adapts to the BH potential, orbital motion, and grid res- 
olution on its (NS) dynamical timescale, i.e., within much less 
than one orbital revolution. This initial phase is almost complete 
by the time the NS crosses the ISCO or becomes tidally dis- 
rupted. 
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The actual values of the mass ratio Q and compactness C 
can be found in Table [1] They are also shown in Fig. |2] and 
were chosen to match closely the models presented in Shibata et 
al. (2009) (see Fig.|3]l. Figures |2] and |3] include the line separat- 
ing the regions within which the NS first reaches the innermost 
stable circular orbit (ISCO) and the mass shedding limit (MS), 
respectively. It is taken from Taniguchi et al. (2008) and is based 
on quasi-equilibrium configurations of the binary system. 

3. Results 

3.1. Resolution and convergence 

Figure |4] shows results for both high-resolution and low- 
resolution models. Simulations with lower mass ratios and less 
compact NSs have intrinsically better resolution (see Table [T]l. 
The dynamics for these models and the results, e.g., gas mass 
on the grid after merger, agree reasonably well for the high- and 
low-resolution calculations. Thus 128^ runs, which are computa- 
tionally feasible in a few days, yield converged results. However, 
the models with very compact NSs - which have comparatively 
fewer zones that cover the NS volume and in particular the steep 
surface gradients - or models with a high mass ratio - i.e., a 
large BH, which in turn implies large zones - require the use of 
at least 256^ zones to be adequately resolved. This can also be 
seen in Fig. |6] for example: for model M2.145, both the low- 
and high-resolution simulations can hold the NS intact on the 
grid and thus the maximum density during the first few orbits 
are equal. On the other hand, the maximum density of the NS in 
model M4.180 is distinctly lower in the low -resolution simula- 
tion than in the high-resolution calculation, which indicates that 
the NS develops a more compact core during the initial orbital 
motion. For very low resolution models, the NS core dissolves 
secularly. For this model the results of the low-resolution simu- 
lation should be interpreted with caution. We have no reason to 
believe that the high-resolution simulations are compromised in 
this way, but are unable to perform even higher-resolution calcu- 
lations to check this point. 

One notices that even for the simulations where the NS-BH 
distance is initially identical, the high resolution models reach 
the point of merger systematically earlier than the low resolu- 
tion cases. This clearly numerical effect seems to decelerate the 
radial inward motion of the lower resolution runs, even in the 
cases resolved best {Q-2) where the central densities for both 
resolutions match (model M2.145 in Fig.|6]l. A detailed analysis 
of the possible cause(s) of this more rapid pre-merging evolu- 
tion of more highly resolved models did not provide an unam- 
biguous conclusion. We could identify differences in neither the 
gravitational-wave energy loss, nor the quality of energy con- 
servation, nor the precision of the gravitational potential energy, 
as clear reason(s) for the observed behaviour. It appears most 
likely that a difference in angular momentum conservation plays 
a role. This quantity, however, is very difficult to calculate pre- 
cisely with its different contributions from gas, orbit, BH spin, 
and gravitational wave losses, and given the model assumptions 
(BH spin effects). Therefore, a quantitatively convincing assess- 
ment cannot be presented for this possibility. 

3.2. Dynamical evolution 

We describe the merger for the two most extreme cases (non- 
spinning BH) in our set of models: (1) a small BH with large 
NS, specifically a mass ratio of Q - 2 and a small value for 
compactness of C = 0.145; (2) a large BH with a small NS, 



2.0 r 




time / ms 

Fig. 6. Maximum density as a function of time for two se- 
lected models. The upper pair of lines shows model M4. 1 80; the 
lower pair of lines shows model M2. 145. Thin lines show lower- 
resolution simulations (128^); bold lines show higher-resolution 
simulations (256-^). 



specifically a mass ratio Q - 5 and a value for compactness of 
C = 0.180. These two models cover the main features seen in the 
simulations, while other models with intermediate values for Q 
and C are easily recognised as being similar to one or the other 
of the two cases. 

3.2.1. Low mass ratio, small compactness 

During the initial orbital decay of model M2.145, the central 
density of the NS remains fairly constant (Figs. |6]and|4|l. One 
notices small oscillations originating in the initially spherical NS 
adapting to the potential of the BH and the discretisation on the 
grid. At about 5 ms, the closest distance of the NS density maxi- 
mum or centre-of-mass is reached (Fig.|5]bottom left; Fig.|7]top 
left). At this point, matter is transferred to the BH, but at a rel- 
atively slow rate compared to the other case. The NS is tidally 
disrupted, but not before a substantial portion of its matter is 
given an elliptic orbit taking it temporarily away from the BH: 
we note the increase in distance of the density maximum in all 
cases of low mass ratio and small compactness (models on the 
lower and left hand side in Fig. |5]l. We return to this point of 
orbital widening in Sect. 14.11 The "radius of the BH" shown in 
Fig.|5]is the arithmetic mean of the event horizon and innermost 
stable circular orbit, as described in Sect. 12.21 For a non-rotating 
BH, this is equal to two Schwarzschild radii. 

The NS matter spreads out in a tidal spiral arc containing 
0.08 Mq (Fig. |7] central left panel r = 9.96 ms). Eventually matter 
falls back along the spiral arc toward the BH. Some of it man- 
ages to complete a full circle around the BH and then collides 
with the incoming material. At this point, the infalling matter is 
deflected away from the BH, the circling material spreads out to 
fill the space and a torus is formed. The point in time when this 
deflection occurs is marked by a short vertical line in Fig.|4] and 
one can see that this applies only to models with low mass ratios 
of Q-2 and Q = 3, for all three values of compactness. In some 
cases, the high-resolution models were not evolved sufficiently 
long to see the formation of a disk. 
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Fig. 4. Temporal evolution of some quantities for various models; the parameter values for Q and C are given within each panel. 
Solid lines show the gas mass. Dotted lines show the maximum density. Thin lines show low-resolution simulations (128^); bold 
lines show high-resolution simulations (256^). A short vertical line indicates the formation of an accretion disk. Short pieces of 
nearly horizontal lines show the amount of unbound mass. 
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Fig. 5. Temporal evolution of some quantities for various models; the parameter values for Q and C are given within each panel. 
The upper pair of lines show binary separation; the lower pair of lines show the radius of the BH (details see text). Thin Unes show 
low-resolution simulations (128^); bold lines show high-resolution simulations (256^). The vertical Une indicates when the mass 
within a sphere around the density maximum falls below 0.03 Mg. 
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A fraction of the material in the spiral arc is unbound. This 
amount is also indicated in Fig. |4] by a short piece of near- 
horizontal line. The values are hsted in Table[T] 

3.2.2. High mass ratio, large compactness 

During the initial orbital decay of model M5.180, the cen- 
tral density increases substantially (Fig. |6] similarly to model 
M4. 1 80) in the calculations of sufficiently high resolution. This 
increase can be noted for all models with massive BHs, i.e., high 
mass ratios (from close inspection of Fig.|4]l. So when the NS is 
close enough to finally transfer mass to the BH, it is very com- 
pact and mass transfer proceeds relatively quickly: we note the 
very steep decline of the mass on the grid for this model M5.180, 
in contrast to M2.145 described above. This effect is also par- 
tially caused by an increase of the mass and radius providing 
an essentially direct approach to the absorbing surface for large 
BHs. 

The distance from the BH at which NS density reaches a 
maximum decreases continuously all the way to the numerical 
surface of the BH. In the right panels of Fig. |5] one can see 
the line of separation that nearly meets the line of BH radius 
to within a few km. The minimal distance shown in the opposite 
case of model M2.145 is much larger, approximately 15 km. 

The NS matter remaining in model M5. 180 after the merger, 
also spreads out in a tidal spiral arc but contains only 10""^ Mq 
(Figs.|4]and[8]l. The BH, which has a very large extent, accretes 
all the material that flows back along the arc; none manages to 
flow all the way around the BH and form a disk. 

Eventually the matter in the outermost portions of the spiral 
arc reaches the edge of the largest grid (only the lower reso- 
lution simulations can be evolved for this length of time). For 
most models, this happens at around 15-20 ms (depending on 
the model) at which point the amount of unbound matter ejected 
from the system is registered by the calculation. The amount of 
matter ejected is approximately 10""'-10"^ Mq. We also calcu- 
late the amount of material on the grid that is not gravitation- 
ally bound, by comparing the internal and kinetic energies with 
the potential. The values of these mass components are listed in 
Table[T]and are consistent with the amounts that actually flow off 
the coarsest grid as described above. Dividing the distance trav- 
elled by this material - roughly half the grid size L/2 «i400 km - 
by the time elapsed between merger and arriving at the edge of 
the grid (approximately 5 ms), yields an average speed of c/4. 

3.2.3. Comparison of models 

Looking at Figs.|4]and|5] one sees that all models with C - 0. 145 
and C = 0.160 follow the patterns described above for model 
M2.145. AU models with C - 0.180 follow the mass accretion 
pattern of M5. 7 SO, although some aspects of models M2.180 and 
M3.180 are different: the BH is small enough in these two cases 
to allow formation of a disk, albeit with a very small mass. We 
included the information about the mass remaining on the grids 
in Fig. |2] in an attempt to outline the similarities between the 
various models. 

3.3. Effects of BH spin 

We calculated four models to assess the influence of BH spin: 
in three models (R2.145, R3.160, R5.180), the BH is taken to 
be initially non-rotating (ai = 0), but we allowed a change of 
the BH potential due to spin-up by accretion of angular momen- 



tum of the NS matter during the evolution. The opposite case, 
an initially rapidly spinning BH with Oi = 0.99, is simulated in 
model S2.145, and like the other three models, the spin and po- 
tential are allowed to change by accretion as well. The complete 
investigation of the effects of these changes is outwith the remit 
of this paper and is left to future work. Here we only point out 
the boundaries of the relevance of results obtained with models 
ignoring BH spin. 

The post-merger rotational speed of the BH increases mono- 
tonically within the three models with initial Oi -Q spin, from 
a value of Of - 0.39 to Of - 0.56 (see Table |2] and Fig. |9]l when 
(decreasing the mass ratio and compactness. Following the pre- 
scription of Eq. ([T]i, an increasing spin reduces both the radius 
of the horizon and the radius of the ISCO (Fig.|9]l. In these mod- 
els, the decrease in radius caused by spin is greater than the in- 
crease caused by mass, or they balance it out, so one observes 
a slight overall decrease in effective radius (model R5.180) or 
hardly any change (models R3.160, R2.145): we note the lower 
near-horizontal line in the right panels of Fig.|9] 

If at the closest point when the NS material is being ac- 
creted by the BH, the radius of the BH remains the same or even 
shrinks — and with it the effective potential — the effective- 
ness of matter accretion is reduced. In the end, after tidal dis- 
ruption, much more material remains in the spiral arc than in the 
equivalent models without BH spin, e.g., compare R5.180 with 
M5.180. When the models that include the effect of BH spin 
are completed, between 0.18 and 0.28 Mq of NS material is 
dispersed around the BH (see Table |2] and Sect. 14.21 for further 
discussion). 

The initially rapidly rotating BH (model S2.145) is deceler- 
ated from a; = 0.99 to Of = 0.88 (Tableland Fig. |9]l because the 
initial BH with ai - 0.99 has a higher specific angular momen- 
tum than the matter that is subsequently accreted. In addition to 
the increase in mass, the reduction in rotation itself increases the 
size of the BH: we note the large increase in the circle radius in 
the right panels in Fig. [8] and the rise of lower solid line in the 
bottom right panel of Fig. |9] 

Model R3.160 is the only one presented in this paper that 
exhibits a distinctive orbit widening, the lower resolution model 
even showing repeated mass accretion episodes. In this case, the 
NS remains self-bound, albeit being less massive, after the first 
close encounter with the BH. In the high-resolution simulation 
the NS is tidally disrupted during the second encounter, while for 
the low-resolution simulation it takes another two episodes be- 
fore the NS is tidally disrupted. The orbital distance plot (Fig.|9] 
right panel, second from top) resembles those from our own 
previous works (Janka et al. 1999) and others (e.g., Davies et 
al. 2005). We return to this point below. 



4. Discussion 

We discuss three specific points: (1) the orbit widening during 
the accretion and tidal disruption of the NS, (2) the amount of 
mass remaining to produce a disk, and (3) the consequences of 
the NS reaching the ISCO before being tidally disrupted. 

4.1. Orbit widening 

In previous simulations that adopted a purely Newtonian poten- 
tial for the BH, it was noted that during mass accretion from 
the NS into the BH, the radial position of the density maximum 
shifted outwards. This resulted in a lighter NS moving away 
from the BH on an elliptic orbit. Of course, the shape of the orbit 
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Fig. 7. The mass density distributions for models M2.145 and M3.180 are displayed in the orbital plane with contours spaced 
logarithmically in steps of 0.5 dex, with units g cm"^. The arrows indicate the velocity field. The circle at the centre outUnes the BH 
radius, which is the arithmetic mean of the event horizon and ISCO. 
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Fig.8. The mass density distributions for models M5.180 and S2.145 are displayed in the orbital plane with contours spaced 
logarithmically in steps of 0.5 dex, with units g cm"^. The arrows indicate the velocity field. The circle at the centre outUnes the BH 
radius, which is the arithmetic mean of the event horizon and ISCO. 
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Fig. 9. Temporal evolu- 
tion for models R5.180, 
R3.160, R2.145, and 
S2.145 (from top to bot- 
tom), which include the 
effect of spin-up of the BH 
(see Table |2]i. The initial 
value of the BH rotation 
parameter is given in the 
panels. Thin lines show 
lower-resolution simula- 
tions (128-'); bold lines 
show higher-resolution 
simulations (256^). The 
left column of panels 
shows: the gas mass (solid 
lines) and the maximum 
density (dotted lines). A 
short vertical line indi- 
cates the formation of 
an accretion disk. Short 
pieces of nearly horizontal 
lines show the amount 
of unbound mass. The 
right column of panels 
shows; the radius of the 
BH (arithmetic mean of 
horizon and ISCO) as 
lower solid line, the binary 
separation as upper solid 
line, and the BH spin 
parameter a as top dashed 
line. The vertical line 
indicates when the mass 
within a sphere around 
the density maximum falls 
below 0.03 Mq. 
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Table 2. Key initial model parameters and some results for models including BH spin. 



Model 


Q 


C 






^NS 


d\ 


A' 
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Ax 


M„ 


Mb 


Mg 


Of 








Mo 




km 


km 




km 


km 


Mo 


Mo 


Mo 


km 


R5.180 


5 


0.180 


7.0 


0.00 


11.5 


88 


256 


1450 


0.71 


0.11 


0.07 


0.18 


0.39 


R3.160 


3 


0.160 


4.2 


0.00 


14.3 


66 


256 


1100 


0.54 


0.08 


0.14 


0.22 


0.44 


R2.145 


2 


0.145 


2.8 


0.00 


14.3 


56 


128 


880 


0.86 


0.05 


0.17 


0.22 


0.56 


S2.145 


2 


0.145 


2.8 


0.99 


14.3 


52 


128 


850 


0.83 


0.02 


0.26 


0.28 


0.88 



Initial parameters: mass ratio Q, compactness C, black hole mass Mbh, initial BH spin parameter a\, neutron star radius i?NS, initial orbital 
distance d\, number of zones per dimension N, size of largest grid L, size of finest zone Ax = L/SN. The mass of the NS is 1.4 Mq in all cases. 
Values at the end of the simulation: unbound+ejected gas mass M^, bound neutron star mass around BH Mb, neutron star mass not instantly 
accreted by BH M^, final BH spin parameter Of. 



and the gravitational wave emission eventually caused the NS to 
experience a second interaction with the BH, in the meantime the 
accretion, however, being effectively shut off. The sequence of 
mass transfer and orbit widening can be repeated. Numerically 
this was observed e.g., in Lee & Kluzniak (1999; their Fig. 20 
shows the binary separation), Janka et al. (1999; their Fig. 1 
shows the orbital separation) and in Rosswog et al. (2004; their 
Fig. 7 shows the mass accretion episodes). It was explained with 
a semi-analytic description by Davies et al. (2005), in which 
some "fraction of the angular momentum of the matter (trans- 
ferred from the NS to the BH) is fed back into the NS". Their 
Fig. 1 shows the orbital separation. 

Miller (2005) argued that the presence of an ISCO, which is a 
purely relativistic effect, would substantially alter the outcome of 
a NS-BH merger, since the NS would be preferentially absorbed 
rather than tidally disrupted. This is basically corroborated by 
in the general relativistic simulations of Rantsiou et al. (2008), 
Etienne et al. (2009) and Shibata et al. (2009; their Fig. 1 shows 
the orbital separation without any sign of orbit widening). 

The Paczyhski-Wiita potential in our simulations mimics the 
ISCO, so the effect oudined by Miller (2005) should be present 
in our otherwise Newtonian simulations. Indeed, the plots of the 
distance between the centre of the BH and the centre-of-mass 
of the NS (Fig. |5] Fig. ^ place the results of our simulations 
squarely in-between the Newtonian and the relativistic results: 
we observe a beginning of an orbit widening, but the widening 
is not fast enough to produce a self-bound 'mini'-NS. In particu- 
lar, in the models with high mass ratios and highly compact NSs, 
the binary distance shrinks to the 'surface' of the BH and hardly 
any subsequent orbit widening is seen. The NS is practically ab- 
sorbed whole, as noted in the relativistic calculations. 

An exception is model R3.160 (including the potential 
changing effects of BH spin), which exhibits the repeated mass 
accretion episodes that are very similar to the purely Newtonian 
models (see Fig.|9]right panel second from top). This might be 
the consequence of two opposing effects. Firstly, the orbit widen- 
ing is more pronounced in models that include BH spin (compare 
right panels of Fig. |9] with the equivalent models in Fig. |5] be- 
cause the denser regions of the 'escaping' NS are caught up to 
a lesser degree within the smaller BH radius. Secondly, model 
R3.160 is intermediate between R5.180 where the large BH 
swallows up the compact NS, and model R2.145 where a small 
BH produces a large differential tidal field on an extended NS. 
In this latter case, the NS is tidally disrupted, although its orbit 
has widened significantly. The remaining mass in model R3.160 
manages to recollect itself into a roughly spherical NS on an 
elliptic orbit. We finally note, however, that the number of mass- 
transfer episodes of this model still varies between the simula- 
tions of different resolution, so the detailed values and sequence 



of events (e.g., number of cycles) have to be interpreted with 
caution. 

4.2. Disk mass 

The amount of material not 'instantly' absorbed by the BH but 
remaining in the surroundings and potentially available to form 
a torus, is an important energy source for powering the jet of a 
gamma-ray burst. Thus the questions of how much matter is ac- 
creted, how much is ejected, and how much forms a torus, need 
to be investigated carefully. In this point, the numerical general 
relativistic simulations (Rantsiou et al. 2008; Etienne et al. 2009; 
Shibata et al. 2009, their Fig. 7) still differ markedly among 
themselves, even for models with similar parameters. The range 
spans from » 0.1 Mq to less than 10^* Mq. To power gamma- 
ray bursts effectively, at least some hundredths of a solar mass 
appear to be needed. 

The values of mass remaining around the BH for our mod- 
els are tabulated in Tables [T] and |2] Only in models with both a 
compact NS (C = 0.180) and a non-rotating BH, does most of 
the matter become accreted quickly by the BH and only 10"^- 
10"^ Mq remain to form a torus, especially for the massive BHs 
with a high mass ratio Q. This trend matches the general rel- 
ativistic findings; however, the absolute quantities of mass are 
higher in our models, e.g., Shibata et al. (2009) obtain typically 
1 % of the neutron star mass for models of low mass ratio Q. On 
the other hand, Etienne et al. (2009) infer 0.06 Mq for a non- 
rotating BH model with Q-3. For less compact NSs, we obtain 
around 0.1 Mq. 

As already mentioned above, our models including the effect 
of BH spin (Table |2]i end with a significant amount (between 
0.18 Mq and 0.28 Mq) of NS material spread out around the 
BH, as the BH's effective potential depth decreases because the 
spin-up matches or outweighs the mass increase (Fig.|9] models 
R5.180, R2.145). 

4.3. Tidal disruption vs. ISCO 

Following Miller's (2005) argument, we would expect little mass 
(tenth of a percent of a solar mass, or less) to remain around the 
BH if the NS reaches the ISCO before being tidally shredded. 
On the other hand, if the NS were tidally disrupted ('mass shed- 
ding', MS) before reaching the ISCO, then significant amounts 
of matter (several percent of a solar mass) would be able to form 
a torus around the BH. With this in mind, and to be able to com- 
pare with previous relativistic simulations, we chose the model 
parameters plotted in Fig.|2] This graph also uses different sym- 
bols to show the varying remaining gas mass around the BH. 
No clear correlation pattern can be seen between the remaining 
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gas mass and the expected line separating the ISCO from the 
MS case, except a general trend that the amount of mass is in- 
deed lower toward the right (models with higher compactness). 
However, we cannot identify a clear step change. This lack of 
clear demarcation between ISCO and MS is also the case for the 
models including general relativity (Shibata et al. 2009, Fig. [3]). 

One point to remember is that the mass ratio Q and compact- 
ness C used in this paper, are based on the purely Newtonian, 
uniquely defined, masses of NS and BH, and radius of the NS. 
Of course, in the models with general relativistic physics, Q and 
C have to be defined with a specific choice of the ADM mass, 
rest mass, circumferential radius, isotropic radius, etc. This com- 
phcates the comparison between our Newtonian results shown 
in Fig.|2]and the general relativistic work of others in Fig. [3] We 
would expect these differences to produce a shift in the separat- 
ing line within the plot, but guess it would only be by several 
tens of percent, and therefore would not qualitatively change the 
statements. 

4.4. Rotating BH 

For models R5.180, R3.160, R2.145, and S2.145, we note only 
that a large difference in dynamics and remaining torus mass is 
present compared to the equivalent non-rotating polytropic mod- 
els (for details see Sect. 13.31 and the end of Sects. l4n and 14. 2l l. 
It is beyond the scope of this paper to map out the differences in 
detail over the full parameter space, but this will be the topic of 
a subsequent investigation. We note that Etienne et al. (2009) re- 
port a fairly high value for the torus mass (0.2 Mq) in their Q-li 
model with a spinning BH (a = 0.75). 

Our values for the final spin parameter of the BH (0.39, 0.43, 
0.56) match the values obtained by Shibata et al. (2009; their 
Table III) fairly well, although their values are systematically 
larger (0.42, 0.56, 0.68, respectively). 

In our pseudo-Newtonian simulations, the mass that assem- 
bles into a torus or is ejected, shows a very strong sensitivity 
to the inclusion of the angular momentum gain by the BH be- 
cause of the accretion of matter in the merger. This effect, how- 
ever, should be generically included in the relativistic models of 
Shibata et al. (2009) and Etienne et al. (2009). 

Since the existence of an ISCO is accounted for by the BH 
potentials used in our simulations, the tendency towards reach- 
ing higher torus and ejecta masses in our models than relativistic 
simulations (with similar system parameters) requires a different 
explanation. We hypothesise that the difference could be a con- 
sequence of the stronger self-gravity of relativistically described 
neutron stars compared to the Newtonian objects considered in 
our simulations. We note how so much more closely the tidal dis- 
ruption limits of the pseudo-TOV models of Ferrari et al. (2009) 
shown in their Fig. 1, match the full GR results than the purely 
Newtonian models. The deeper relativistic potential might im- 
pede the disruption of the neutron star and the formation of a 
tidal arc, so that a greater fraction of the star is directly swal- 
lowed by the BH. 

5. Conclusions 

We reach the following conclusions. 

1. We do not note any fundamental and discontinuous differ- 
ence in the dynamics of the models in the ISCO regime as 
opposed to the 'mass shedding' (= tidal disruption) regime. 
The variation in masses that are not instantly accreted by the 
BH is continuous (Figs.|2]and[3]i. 



2. In all cases the NS is initially tidally stretched into an arc, 
albeit of very different mass in each of the various cases. 

3. For models without BH spin, the onset of orbital widening 
is noted during the phase of tidal shredding of the NS but is 
only short-lived. 

4. In one case of a model including the BH spin, we see the for- 
mation of a 'mini' NS on an extended elliptic orbit, followed 
by episodic mass transfer, and a final tidal disruption. 

5. In mergers with compact neutron stars (C = 0.180), a very 
small amount, less than lO^^M©, remains in the surroundings 
of the BH. 

6. Only in mergers of low mass BHs with mass ratio Q-2 or: 
Q-J>, does a closed accretion torus - as opposed to an open 
arc - form. 

7. Mergers with neutron stars of compactness C - 0.145 or 
C - 0.160 produce a significant amount of material (order 
of 0.1 M0) in the surroundings of the BH. 

8. Effects due to BH spin-up are important and significantly 
change the results. We note an increase in mass of the mate- 
rial surrounding the BH of up to 0.2-0.3 Mq. 
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